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Abstract 

We identify an intriguing feature of the electron- vibrational dynamics of molecular systems via 



a computational examination of trans-p olyacetylene oligomers. Here, via the vibronic interactions, 



£> \ the decay of an electron in the conduction band resonantly excites an electron in the valence 

00 

CO ' band, and vice versa, leading to oscillatory exchange of electronic population between two distinct 



electronic states that lives for up to tens of picoseconds. The oscillatory structure is reminiscent 

o 

CO 



£^ , of beating patterns between quantum states and is strongly suggestive of the presence of long- 



lived molecular electronic coherence. Significantly, however, a detailed analysis of the electronic 
coherence properties shows that the oscillatory structure arises from a purely incoherent process. 
These results were obtained by propagating the coupled dynamics of electronic and vibrational 
degrees of freedom in a mixed quantum-classical study of the Su-Schrieffer-Heeger Hamiltonian 
for polyacetylene. The incoherent process is shown to occur between degenerate electronic states 
with distinct electronic configurations that are indirectly coupled via a third auxiliary state by the 
vibronic interactions. A discussion of how to construct electronic superposition states in molecules 
that are truly robust to decoherence is also presented. 
I. Franco, A. Rubio and P. Brumer New J. Phys. 15, 043004 (2013) 
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I. INTRODUCTION 



A characteristic feature of molecular systems is that they exhibit strong electron- 
vibrational interactions. Such vibronic couplings |l|, |2[ are an essential component of the 
photophysics of molecules, leading to vibrations upon electronic excitation [3[ , spectral line 
broadenings, nonradiative transitions |4j-[7|, electronic relaxation [8[ and decoherence [9HL4]. 

In this paper we identify an intriguing and novel feature of the electron-vibrational dy- 
namics of trans -polyacetylene in which, via the vibronic interactions, the decay of an electron 
in the conduction band leads to resonant excitation of an electron in the valence band. The 
converse process (the decay of an electron in the valence band to a further inner state leading 
to excitation of an electron in the conduction band) also takes place and brings the system 
back to its original state. The result is long-lived oscillatory electron dynamics. Through- 
out we refer to this phenomenon as Vibronically-Induced Resonant Electronic Population 
Transfer (VIBRET). 

As a model of trans-polyacetylene (PA) we employ the Su-Schrieffer-Heeger (SSH) Hamil- 
tonian 15] , a tight-binding model for PA with strong electron-vibrational interactions. The 
SSH Hamiltonian is often used to study the static and dynamic features introduced by 
strong electron-ion couplings in molecular systems 



12 



Il6l . Il7| . It has been shown to be 



successful in capturing the basic electronic structure of PA, its photoinduced vibronic dy- 



namics and the rich photophysics of polarons, breathers and kinks [3|, ll2|, ll8|, |l9|. The 



vibronic dynamics of SSH chains is followed by explicitly propagating the coupled dynamics 
of electronic and vibrational degrees of freedom in an Ehrenfest mixed quantum-classical ap- 



proximation 12 



20M22| where the nuclei are treated classically and the electrons quantum 



mechanically. Effects of nuclear fluctuations and decoherence are captured by propagating 
an ensemble of trajectories with initial conditions obtained by sampling the ground-state 
nuclear Wigner phase-space distribution J23]. 

Below we show that for a specific class of initial states this model of the vibronic evo- 
lution leads to VIBRET in SSH chains that, depending on system size, can live for up to 
tens of picoseconds. VIBRET is seen to arise between degenerate electronic states with 
distinct electronic configurations that are indirectly coupled via a third auxiliary state by 
the electron-vibration interactions in the system. Given this identified level structure, we 
investigate the effect of changing system size and the nature of the initial state on the 



dynamics. 

A striking feature of the VIBRET is that it leads to population oscillations among the 
relevant levels that are analogous to those observed in beatings that result from coherent 
superposition states. As such, these oscillations seem to indicate that underlying this dy- 
namics is an electronic superposition state that can live for picoseconds, a timescale that is 

24| . The question of whether the observed be- 



very long for electronic coherences 
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havior is, in fact, due to a long-lived electronic coherence is particularly relevant because of 
spectroscopic observations in photosynthetic systems that suggest that unusually long-lived 



electronic coherences are possible in the Fenna-Matthew-Olson and related complexes 25 



271 ]. with timescales exceeding 400-600 fs. Such long-lived electronic coherences have also 



been noted in intrachain energy migration in conjugated polymers 28J. Hence, if the SSH 
model can sustain long-lived coherences even in the presence of strong vibronic couplings, 
an analysis of the coherence properties of the model may well shed light on this topical 



problem 29 



-m- 

We have divided this analysis into three main components. In Sec. IIII Al we discuss the 
essential phenomenology of VIBRET and clarify the basic structure behind the population 
oscillations. In Sec. IIII Bl we characterize the coherence properties of VIBRET by introducing 
reduced measures of the purity that are apt for many-particle systems. Using these purity 
measures we show that, contrary to intuition, the long-lived oscillations observed during 
VIBRET are the result of an incoherent process. Last, in Sec. IIII CI we demonstrate how to 
construct electronic superpositions in vibronic systems that are truly robust to decoherence. 
This set of results is expected to have implications in our understanding of vibronic and 
coherent-like phenomena in molecules, macromolecules and bulk materials. 

II. MODEL AND METHODS 

In the SSH model, the PA is described as a tight-binding chain, where each site represents 
a CH unit, in which the 7r-electrons are coupled to distortions in the oligomer backbone by 
a parametrized electron-vibrational interaction. For an A-membered oligomer, the SSH 



Hamiltonian reads 15 ] 



Hssh = H e + H pix , (la) 



where 

AT-l 
He = ^ i~ f + a ( U n+l ~ Un)]{ci +ltCT C n>ff + cl i(T C n+l!(T ) (lb) 

n=l,ff±l 



AT N-l 



^ = E^ + f E(^- m «) 2 ' ( lc ) 

n=l n=l 

are, respectively, the electronic (H e ) and nuclear (H p h) parts of the Hamiltonian. Here 
u n denotes the displacement of the nth CH site from the perfectly periodic position x = 
na with a the lattice constant of the chain. M is the mass of the CH group, p n is the 
momentum conjugate to u n and K is an effective spring constant. The operator c [ na (or 
c„ iCr ) creates (or annihilates) a fermion on site n with spin o and satisfies the usual fermionic 
anticommutation relations. The electronic component of the Hamiltonian consists of a 
term describing the hopping of 7i electrons along the chain with hopping integral to an d 
an electron- ion interaction term with coupling constant a. The quantity a couples the 
electronic states to the molecular geometry and constitutes a first-order correction to the 
hopping integral that depends on the nuclear geometry. Throughout this work, we assume 



neutral chains with clamped ends and use the standard set of SSH parameters for PA [151 ]: 
t = 2.5 eV, a = 4.1 eV/A, K = 21 eV/A 2 , M = 1349.14 eV fs 2 /A 2 , and a = 1.22 A. While 
it is possible to supplement the model with on-site electron-electron interaction terms, for 
the discussion below these terms are not fundamental and do not change the main findings. 
We therefore focus on the usual case of noninteracting electrons coupled to phonons. 

The method employed to propagate the electron-vibrational dynamics of SSH chains 
has been described in details previously [12|, [14|. Briefly, the dynamics is followed in the 



Ehrenfest approximation 20], where the nuclei move classically on a mean- field potential 



energy surface with forces given by 



, 9Hs$ B 



Pn = ~(<P(t)\-^Ht)). (2) 

In turn, the antisymmetrized many-electron wavefunction \<p(t)) satisfies the time- dependent 
Schrodinger equation 

ih^\<p(t)) = HMt)Mt))> (3) 

where u = (ui, U2, ■ ■ ■ , %). Decoherence effects are incorporated by propagating an ensem- 



ble of quantum-classical trajectories with initial conditions se 
of the ground-state nuclear Wigner distribution function 



ected by importance sampling 
23J of the oligomer obtained 



in the Harmonic approximation. In this way the dynamics reflects the initial nuclear quan- 
tum distribution and is subject to the level broadening and internal relaxation mechanism 
induced by the vibronic couplings. Results shown here correspond to averages over 10000 
trajectories, providing statistically converged results. 

In its minimum energy conformation, the SSH Hamiltonian yields a chain with a per- 
fect alternation of double and single bonds. Its electronic structure is composed of N/2 
"valence band" orbitals with negative energies and N/2 "conduction band" orbitals with 
positive energies that in the long-chain limit are separated by an energy gap of 1.3 eV. The 
single-particle spectrum depends on the nuclear geometry and changes during the electron- 
vibrational dynamics. However, because of the electron-hole symmetry in the Hamiltonian, 
the orbital energies are always such that for each orbital in the valence band of energy — q 
there is an orbital in the conduction band of energy 6j. 

III. RESULTS AND DISCUSSION 
A. VIBRET 

1. Basic phenomenology 

We begin by describing the basic phenomenology behind VIBRET. For this, consider the 
dynamics of an oligomer initially prepared in a separable superposition state of the form 

\tt) = {b \%) + b 1 \<$> 1 ))®\xo) (4) 

where the |xo) is the initial nuclear state that, for definitiveness, we take it to be the ground 
nuclear state of the ground electronic surface. As a first example, consider bo = b\ and |$o) 
and |$i) to be the states obtained by HOMO-^LUMO and HOMO^LUMO+1 transitions 
from the ground state in a given spin channel. That is, 

l$o) = Clumq^chomcvIvo) 

t (5 ^ 

where \tpo) is the ground electronic state. Because of the electron- hole symmetry in the SSH 
approach, state |$i) is degenerate with the state 

l $ 2> = c[ UMO C H OMO-l,a|y?0>- (6) 
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FIG. 1: VIBRET in SSH chains with 4, 20 and 100 sites initially prepared in the superposition 
defined by Eqs. (jl])-©, with bo = b±. The figure shows the population of the orbitals involved 
during the complex vibronic evolution. Note the electronic population exchange among levels for 
N = 4 and N = 20. 



States |<&i) and |$2) have distinct electronic configuration and are thus orthonormal. 

Figure Q] shows the population dynamics in the four relevant orbitals (HOMO-1, HOMO, 
LUMO, LUMO+1) involved during the dynamics for chains of different size (N = 4, 20 
and 100) . Focus first oniV = 4 (Fig. [T]A.) where the vibronic dynamics shows a remarkable 
long-lived population transfer between levels that survive for tens of picoseconds. Such 
population exchange is what we refer to as VIBRET. An analysis of the orbital dynamics 
leads to the following interpretation of the observed behavior (a schematic diagram of the 
population exchange in a generalized setting, and from a single-particle perspective, is shown 
in Fig. [2]). The population from the higher energy conduction band orbital f (in Fig. QJ 
f =LUMO+l) is transferred into the lower energy conduction band orbital i' (i' =LUMO in 
Fig. [1]). Through the vibronic interactions, this decay resonantly drives an electron from the 
inner-most valence orbital j (j =HOMO-l in Fig. [1]) into the higher energy valence orbital i 
(i =HOMO in Fig. CD). Since the orbital energies are such that ey — ey = e« — e 7 , the process 
is energy conserving. In the case shown in Fig. [TK, complete population transfer occurs in 
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FIG. 2: Scheme of the single-particle states and the dynamics involved during the VIBRET. In the 
figure the horizontal lines denote orbital levels of varying energy E and the circles denote electrons. 
The arrows indicate the joint electron exchange observed during VIBRET: the decay of an electron 
from the higher-energy conduction band orbital j' into a lower energy state i! leads to resonant 
excitation of an electron from an inner state in the valence band j to a higher energy state in the 
valence band i. Upon population inversion, the electron in the valence band i decays back into 
state j and resonantly excites the electron in level i! to level j' . Several of these cycles can be 
observed when the population exchange is energy conserving, i.e. when ej — €j = e' — e[. The 
complementary iV-particle level structure and couplings are depicted in Fig. El 
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FIG. 3: Level structure and couplings between iV-particle electronic states able to sustain VIBRET. 
The vibronic dynamics couples the degenerate states |3>i) and \&2) indirectly through a third 
auxiliary state |$o) via vibronic nonadiabatic coupling terms Vna [see Eq. (jHJ)] forming a A or V 
level system. Figure [2] depicts one possible single-particle electronic distribution for the states in 
the triad. 



~3 ps. At this stage, as shown in Fig. |2j the converse process occurs. Transfer of population 
from state i to the lower energy state j in the valence band resonantly excites population 
from i' into j' in the conduction band. From the perspective of many-particle states, the 
system is effectively transferring population from state |$i) to state ($2); and vice versa. 
The computed vibronic evolution for the four site chain (Fig. [TJA.) reveals at least six of these 
|$i) — > I $2) — > l^i) cycles before the population equilibrates. 

We have also observed this intriguing oscillatory population dynamics in 20-site chains 
(Fig. [IB)- Here, the period of oscillation (~ 0.7 ps) is faster than in the N = 4 case and 
observable for 2 to 3 ps. That is, while the population dynamics is evident for N = 20, the 
process competes with population transfer into other electronic states that are coupled by the 
vibronic evolution. Such competition substantially decreases the lifetime of this population 
exchange with respect to the N = 4 example, where levels are spectrally isolated. In fact, 
for iV = 100 (Fig. [Tp) the electronic spectrum is so dense that the VIBRET is simply not 
observed. 

2. Level structure and underlying couplings 

A closer look into the three iV-particle eigenstates of H e involved in the dynamics, |3>o)> 
|$i) and I $2)) reveals that they conform to the three-level system schematically shown in 
Fig. El Such system consists of two degenerate orthonormal states, |$i) and ($2), with 
distinct electronic structure that are coupled indirectly through a third eigenstate |$ ) via 
nonadiabatic coupling terms Vna- The states |$i) and |$ 2 ) are uncoupled in the vibronic 
evolution and thus require of the third "auxiliary" state |$o) i n order to transfer population 
between one another during the electron-vibrational dynamics. This resonant population 
transfer is a second order process in Vna that leaves the population of the auxiliary state 
approximately constant through the dynamics. 

In order to understand how these effective couplings between levels |$o) 5 l^i) an d ($2) 
arise, consider the selection rules for nonadiabatic couplings in the context of (generic) mixed 
quantum-classical dynamics. For an electronic wavefunction \<p(t)) that satisfies Eq. ([3]) 
where the time-dependence in H e (u) is assumed to arise from the fact that the nuclei satisfy 
some trajectory u(t), the coefficients in the expansion of \f{t)) in terms of adiabatic eigen- 



states, i.e. \ip(t)) = J2 k c k{t)\v k [u(t)]) where H e (u)\<p k [u(t)]) = E k (t)\cp k [u(t)]), satisfy |43j]: 



ih—^ = Ei(t)ci - ih^2(ipi\ — \ip k )c k . (7) 

jfc 

The second term is the nonadiabatic coupling V ik between adiabatic states i and k and can 

be expressed as: 

Vifc = -inu ■ {Pi\\/ U \(p k ) = ihu . (8) 

m — &k 

Here we have taken into account the fact that the sole time dependence of the adiabatic 
states is through the nuclear coordinates such that §i\tp k ) = u ■ V u \(p k ). In SSH chains, 
V ik 7^ only if (ipi\V u H e \ip k ) ^ 0. This is the case, even for degenerate states. Since H e is a 
single particle operator, it then follows that Vi k ^ between states that differ by at most a 
single-particle transition. Two particle transitions require terms in the Hamiltonian that are 
quartic in the creation and annihilation operators, that are absent from this model. More 
generally, for systems in which the electron-phonon coupling term in the Hamiltonian is 
quadratic in the creation and annihilation operators, in order for states to conform to Fig. [3] 
it suffices to guarantee that the selected states |$i) and |<&2) are degenerate and differ by 
two particle transitions, but that they are both a single-particle transition away from some 
state |$o)- This holds even in the presence of electron-electron interactions as they do not 
contribute to the nonadiabatic transitions. The states employed in Fig. Q] satisfy precisely 
these requirements. 

It is now natural to ask if other triads that conform to the scheme in Fig. |3]will also display 
VIBRET. For illustrative purposes, we focus on the N = 4 example. In this case there are 19 
possible electronic states and 5 possible degenerate manifolds (without taking into account 
spin degeneracies); they are tabulated in Table [H The states are labelled by the population 
of its four eigenorbitals, in ascending order. In this notation, the ground state would be state 
(2200), first excited state (2110), etc. For instance, the superposition employed in Fig. [QA. 
would correspond to (2101) + (2110). Figure H] shows the orbital populations during the 
dynamics of chains initially prepared in the state of Eq. (j4]) with 6o = b\ and for different 
choices of |$o) and |$i). The two states involved in the superposition are indicated in each 
panel. The auxiliary state |$o) is further labeled by an 'a' after the orbital occupations. 
Figure H]A. corresponds to a situation similar to that described in Fig. [TJ\ and Eq. (jSJ) . The 
auxiliary state |$ ) is the same but now the population, instead of being initially in state 
(2101), is initially allocated to the second state in the degenerate manifold (1210). A similar 



dynamics results, confirming the basic identified level structure. The dynamics exemplified 
by Fig. [4j3 also uses the (2101)-(1210) degenerate manifold but here the auxiliary state 
|$o) is higher in energy, forming a A system instead of the V system explored in Fig. [TJA. 
Long-lived population transfer between |$ x ) and |$2) is also evident in this case but with a 
different timescale resulting from the change in the nonadiabatic coupling due to change in 
|$o)- Figures HP and D demonstrate the effect in higher energy degenerate manifolds. In 
all cases considered population transfer is as described in Figs. [2] and [31 and survives for 
tens of picoseconds. Naturally, if the system is prepared in a state that does not conform 
to the scheme in Fig. [3] no VIBRET will result. Figure SI exemplifies such a situation for a 
system initially prepared in state (2200) + (2110). 

Note that because the population exchange occurs between degenerate electronic states, 
there is no net absorption or emission of real phonons during the process. Hence, the 
limitations of Ehrenfest dynamics in describing spontaneous emission of phonons 16J do not 
play a significant role here. 
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TABLE I: Electronic distribution and energy of all 19 possible states for a N = 4 BA chain. 
The distribution correspond to the occupation of the four eigenorbitals in ascending energy. The 
electronic energy is given at the ground state optimal geometry {u\ = U4 = 0, U2 = 0.0847 A, 
us = — m). Note the degeneracies in the electronic spectra. 
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FIG. 4: Different types of initial electronic superpositions |0) = (l/v / 2)(| < l ) o) + l^i)) ® Ixo) f° r 
N = 4 that exhibit the VIBRET. The distribution of the states involved in each case are: (A) 
|*o) = (2110), |$i> = (1210); (B) |* ) = (1201), |$i) = (2101); (C) |$ ) = (2020), |$i) = (2011); 
(D) |$o) = (0202), |$i) = (1102). The energies of such states are shown in Table HI The numbers 
in the legend correspond to the eigenorbital labels in ascending energy. 

3. Dependence on the amplitudes of the initial superposition 

Another significant aspect of the observed behavior is the dependence of the VIBRET 
on the amplitudes of the states involved in the superposition at the time of preparation. 
Figure shows the orbital population dynamics for a chain with 4 sites initially prepared 
in the superposition defined by the states in Eqs. ©-© for different bo and b±. When all 
the population is in the auxiliary state (i.e., |6 | 2 = 1) n o population transfer is observed 
(see Fig. [5JA) because the vibrational degrees of freedom are not able to resonantly couple 
the auxiliary state |$o) with the degenerate manifold. As the population initially placed 
in the excited state manifold is increased (the progression shown in Fig. |5B-F) the amount 
of population exchanged during the dynamics changes. Because of the resonance structure 
in Fig. [31 only the population that is initially placed in the degenerate manifold can be 
exchanged, e.g. the state in which initially |&i| 2 = 0.9 can exchange at most 0.9 electrons. 
In addition, by changing the initial coefficients in the superposition one is changing the forces 
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FIG. 5: Dependence of the VIBRET on the amplitudes of the initial superposition. In this example, 



the initial state is given by |0) = (&o|$o) +&l|^i))IXo) with bo = \/l — |&i| 2 . Results are for N = 4 
and for the superposition defined by the states in Eq. ([5]), i.e. |$o) = (2110) and |<I>i) = (2101). 



They correspond to (A) |6i| 2 = 0; (B) |&i| 2 = 0.1; (C) |&i| 2 = 0.25; (D) |&i| 2 = 0.75; (E) \bi 



0.9; 



(F) |6i| 2 = 1. The case of |6o| 2 = |&i| 2 = 0.5 is shown in Fig. Q]A. In all cases the initial coefficients 
where chosen to be real and positive. The numbers in the legend correspond to the eigenorbital 
labels in ascending energy. 

that act on the nuclei at t — and thus, effectively, the strength of the nonadiabatic coupling 
terms during the evolution. This change in the strength of the nonadiabatic couplings 
modifies the timescale of the population oscillations and the lifetime of the process. Note 
that when all population is placed in the degenerate manifold (Fig. [5]F) the strength of 
the nonadiabatic coupling terms is substantially reduced and no population exchange is 
observed in the simulated time window. We have observed this behavior in all of the cases 
considered and, as such, it is an inherent feature of this highly nonlinear vibronic evolution. 
In principle, the phenomenon does not require an initial superposition between states |$o) 
and |$i). In practice, however, such initial coherences introduce additional forces on the 
nuclei at initial time that enhance the effective non-adiabatic couplings, leading to a visible 
effect within the propagated time window. 



12 



B. Electronic coherence during the VIBRET 

Consider now the electronic coherence properties of the VIBRET in order to determine 
whether the observed dynamics is coherent or incoherent. For definitiveness we focus on 
the dynamics in Fig. [TJA. in which the system is initially prepared in the state defined by 
Eqs. (JH)-(jSD with N = 4 and bo = b\. The top panel of Fig. [6] shows the populations of 
the states |$o)j |$i) an d l^), reconstructed from the orbital populations by supposing that 
only these three many-particle states participate in the dynamics. The plot clearly shows 
the population exchange between |$i) and ($2), and its decay, while the population of state 
|$o) remains approximately constant throughout. 

The dynamics exemplified in the top panel of Fig. [H] (also Fig. [TJA.) strongly suggests the 
presence of a long-lived electronic coherence because they are reminiscent of beating patterns 
resulting from superpositions between nearly degenerate states. If, in fact, the dynamics is 
a coherent process then the observed evolution would constitute a clear example of a long- 
lived coherence that is unquestionably electronic. That is, here, the observed beatings 
could arise from the effective coupling between states |$i) and [$2) that is introduced by 
the nonadiabatic coupling terms. The decay of the population exchange in the degenerate 
manifold would then suggest that a decoherence process is taking place with an unusually 
long decay constant (typical decoherence timescales obtained with this model of the vibronic 
evolution are of ~10-100 fs [141]). 

To examine this possibility we now quantify the coherence properties during VIBRET. 
To proceed, it is useful to recall some basic facts about electronic decoherence in molecular 
systems. Electronic decoherence in molecules arises because of interactions with the nuclear 
degrees of freedom and can be understood in terms of nuclear dynamics on alternative elec- 
tronic potential energy surfaces [111 . 1 121 . uA |44j . To see this, consider the reduced electronic 



density matrix associated with a general entangled vibronic Born-Oppenheimer state of the 
form \n(t)) = J2 n e- tEnt/h M\Xn(t)), 

p e (t) = Tr N {\Q(t))(Q(t)\} 

nm 

Here the trace is over the nuclear states, \(p n ) are the electronic eigenstates [H e i ec \(p n ) = 
E n \<f n )], \xn{t)) the nuclear wavepacket associated with each electronic level and u nm = 

13 



(E n — E m )/h. Note that the magnitude of the off-diagonal elements of p e (t) are proportional 
to the nuclear overlaps S nm (t) = (Xm{t)\Xn{t)) ■ Hence, the loss of such p e {t) coherences is 
a result of the evolution of the S nm (t) due to the vibronic dynamics. Standard measures of 
decoherence capture precisely this. For example, the purity of such entangled vibronic state 



45l | is given by 

Tr(^(t)) = ^|( Xm (t)|Xn(t))| 2 (10) 

urn 

and decays with the overlaps of the nuclear wavepackets in the different electronic surfaces. 
In order to quantify the coherences during VIBRET one ideally would like to study the 
purity [Eq. (TTOj) ] directly. However, for many-electron systems, like the one considered here, 
the electronic density matrix p e is a many-body quantity that is not easy to compute and 
hence reduced descriptions of the purity are required. Here we introduce and follow the 
dynamics of the one-body and two-body reduced purities, defined as: 

P 1 (t) = Tr{p 2 (t)}, 
P 2 (t) = Tr{f 2 (t)}, 

where p and T refer to the one-body and two-body electronic density matrices. These 
quantities are defined as: 

p£ = ^Tr{cJ ff cv/) e } (12) 

a 

T pq = 2 Yl Tt { c IA<t ,c ™' Cs °m- ( 13 ) 

a.a 1 

Because the one-body purity is constructed from the one-body density matrix, it only in- 
forms about coherences between states that differ at most by single-particle transitions. For 
example, it cannot distinguish between a superposition and a mixture between states that 
differ by two (or more) particle transitions. Similarly, the two-body purity is only informa- 
tive about the coherences between states that differ by at most two particles transitions. 

The middle and bottom panels of Fig. [6] shows the dynamics of the reduced purities 
during VIBRET. In order to interpret the results we consider three models of the state of 
the system that differ in the assumed degree of coherence. The models are defined by the 
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FIG. 6: Electronic decoherence during VIBRET. The plot shows the populations of the many 
body states |$o)> l^i) and |3>2) (top panel), the one-body purity (middle panel) and the two-body 
purity (bottom panel) for a N = 4 system initially prepared in the state defined by Eqs. dH)-© 
with 60 = b\. The insets highlight the first 200 fs of evolution which are not resolved in the main 
plots. In the purity plots, the simulated data is shown in black. The colored lines represent a fully 
coherent (Ml), partially coherent (M2) and completely incoherent (M3) model of the state of the 
system during the dynamics. In the middle panel the M2 and M3 lines are on top of one another 
and cannot be distinguished. 

following assumed forms for the electronic density matrix: 



Ml: p e = (c |$o) + Ci|$i) + c 2 |$ 2 ))(c*($ | + c*($i| + c*($ 2 |); 
M2: p e = |c | 2 |$ )($o| + (ci|$i) + c 2 |$ 2 ))(ct($i| + c*($ 2 |); 
M3: Pe = |c | 2 |$ )($o| + |ci| 2 |$i)($i| + |c 2 | 2 |$ 2 )($ 2 |. 



(14) 



They represent, respectively, a fully coherent model (Ml), a partially coherent model where 
only the coherences within the degenerate manifold are maintained (M2), and a fully inco- 
herent model (M3). The purity resulting from the simulation is shown in black, while the 
colored lines show the purity expected for these three models, reconstructed by supposing 
that only the |3>o)? I^i) an d |$ 2 ) many-particle states participate in the dynamics. The 
insets highlight the first 200 fs of evolution which are not resolved in the main plots. 

From Pi(t) (middle panel) we conclude that a fully coherent picture is not representative 
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of the actual dynamics. Here, the system begins in a pure state and during the first 200 fs 
the system displays fast decoherence between the states in the initial |<&o)>|$i) superposition. 
The inset details this initial decoherence process. The recurrences observed in the one-body 
purity signal the vibrational dynamics in the excited state manifold [14j that lead to time- 
dependence of the overlaps in Eq. (1101) . These recurrences are not captured by the models 
in Eq. (Tl4|) because they do not take into account the nuclear evolution. After this initial 
fast decoherence, P\ oscillates, reflecting the population changes in the system throughout 
the dynamics. However, Pi cannot distinguish between the partially coherent model M2 and 
the fully mixed case M3 because the coherence in M2 is between states that differ by two- 
particle transitions. In order to distinguish between these two cases we follow the two-body 
purity P2{t) shown in the bottom panel. This quantity shows an initial fast decay (in ~ 
200 fs) due to the decoherence between states |$o) and |$i), followed by oscillations. After 
this initial fast decoherence dynamics the model that best adjusts to the observed behavior 
is M3. That is, the observed population exchange, even when reminiscent of beatings in 
coherent superpositions, is really best described as a mixed state between |$o) 5 l^i) and 
| $2)- We thus are forced to conclude that, contrary to intuition, after 200 fs the dynamics 
during VIBRET is a purely incoherent process. 

The dynamics is incoherent because superpositions between states |$o) and |$i) and |$o) 
and I $2) decohere quickly (on the order of ~200 fs in the N = 4 case and of 10s of fs for 
larger oligomers). Since all communication between |$i) and |$2) is then through |$o) 5 a 
net incoherent process results. From a quantum-classical perspective, the Vna in individual 
trajectories leads to coherences between the states. However, in average the Vna lead to an 
incoherent coupling contributing to the incoherent dynamics. The observed directionality in 
the population exchange is due to the population imbalance between states |$i) and |$2)- 
In fact, had we started with a state where |$i) and |<J> 2 ) where equally populated, then no 
VIBRET would result. 
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FIG. 7: Two-body purity for aJV = 4 chain prepared in state Eq. (|15p . The quantity P2W shows 
that this state starts and remains pure during the dynamics even when it is not a stationary state 
of the Hamiltonian. 

C. What would constitute a superposition of electronic states that is truly robust 
to decoherence due to vibronic couplings? 

Consider the dynamics where the same vibrational wave packet is prepared in two degen- 
erate electronic states, that is 

|n> = -^(|*i> + |* 2 »®lx>- (is) 

Figure [7] shows the time-dependence of the two-body purity for a N = 4 system prepared in 
Eq. ( IT5l) with |$i) and ($2) defined by Eqs. (jSD and (jEJ). For a coherent superposition one 
expects Pi(f) ~ 5.0 in this case, while a perfectly incoherent state would yield Pi{t) = 4.5. As 
can be seen, even when the initial state in Eq. ( 1T5|) is not an equilibrium state of the vibronic 
Hamiltonian and leads to a complex electron-vibrational evolution, this initial superposition 
starts and remains pure throughout the dynamics. It constitutes a clear example of an 
electronic superposition state with coherence properties that are robust to the vibronic 
interactions of the chain. 

The feature that underlies these robust coherences is the fact that the superposition is 
between two electronic states with underlying potential energy surfaces that differ at most by 
a constant factor. In this case, the states |$i) and |$2) are degenerate in all conformational 
space, i.e. E\{u) = ^(w) for all u, where E{(u) is the potential energy surface associated 
with adiabatic state i. Consequently, a given vibrational wavepacket |%) will move identically 
on both surfaces and the nuclear overlap S\2(t) = (xi(t)\X2(t)) that determines the electronic 
coherence between the two states [recall Eq. (jUJ)] is unaffected by the dynamics. In other 
words, the two levels involved couple to the environmental bath identically, preventing the 
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bath from entangling with (and thus inducing decoherence of) the system. This is the case 
provided that the two states are spectrally isolated from other electronic states. 

The quantum structure involved in these robust coherences falls into the class of deco- 
herence free subspaces [46[ |47|. Such subspace has been suggested (see, e.g., Ref. [48!]) to 
underlie the long coherences in the photosynthetic example, although approximate explicit 



computations 



49| have not yet revealed such a structure. By contrast, the triad in Fig. [3] 
does not conform to a decoherence free subspace because the potential energy surface of the 
third electronic state |<J>o) generally differs by more than a constant to that of the degenerate 
states. 

IV. CONCLUSIONS 

We have identified a new basic feature of the vibronic evolution of a molecular system 
that we term VIBRET (Vibronically Induced Resonant Electronic Population Transfer). 
In this process, via the vibronic interactions, the decay of an electron in the conduction 
band to a lower energy state resonantly excites an electron in the valence band, and vice 
versa. In PA oligomers (as described by the SSH Hamiltonian in a mixed quantum-classical 
approximation) the population transfer can survive for up to tens of picoseconds and ob- 
serve several cycles of population exchange. The process requires two degenerate electronic 
states with distinct electronic configurations that are indirectly coupled to a third state via 
vibronic interactions. For Hamiltonians with electron-phonon coupling terms that are at 
most quadratic in the fermionic operators, such population exchange is realized between 
degenerate states that differ by two particle transitions but that are both a single-particle 
transition away from a third auxiliary state. 

The observed population dynamics is strongly suggestive of an electronic coherent process 
with an unusually long decoherence time. However, things are not always what they seem 
and, contrary to intuition, an analysis of the one-body and two-body electronic purities 
shows that VIBRET occurs incoherently. 

We have also demonstrated electronic superpositions in a molecular system that is robust 
to decoherence induced by vibronic couplings. As shown, robust electronic superpositions 
can arise when the underlying potential energy surfaces of the states involved in the su- 
perposition differ by a constant factor. Under such conditions the vibronic evolution of an 



18 



initially separable state does not lead to entanglement between the electronic and vibrational 
degrees of freedom and thus does not lead to decoherence. 

We expect the phenomena described here to be of importance in the understanding of 
vibronic and coherence phenomena in molecules, macromolecules and bulk materials. Future 
prospects include performing fully quantum simulations of VIBRET and determining ways 
to manipulate the identified robust electronic coherences. 
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Fig. S 1: Orbital population dynamics for a 4 site SSH chain initially prepared in a superposition 
with equal coefficients between the ground (2200) and first excited (2110) electronic state. The 
initial nuclear state is taken to be the ground vibrational state of the ground electronic surface. 
This superposition does not conform to the scheme in Fig. 3 and hence does not exhibit VIBRET, 
as evidenced by the constant orbital populations throughout the vibronic evolution. The numbers 
in the legend correspond to the eigenorbital labels in ascending energy. 
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